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Abstract 

We argue that collinearly ordered states which exist in strongly frustrated spin 
systems for special rational values of the magnetization are stabilized by thermal as 
well as quantum fluctuations. These general predictions are tested by Monte Carlo 
simulations for the classical and Lanczos diagonalization for the S = ^ frustrated 
square-lattice antiferromagnet. 
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Frustrated 2D quantum spin systems can give rise to interesting behavior in 
high magnetic fields. For example, plateaux in the magnetization curve of 
SrCu 2 (B0 3 ) 2 have been observed at0 (M) = §, 7 and | [1]. This has sparked 
theoretical interest in the magnetization process of the related 2D Shastry- 
Sutherland model (see e.g. [2]). 

Owing to its simplicity, the frustrated square-lattice antiferromagnet (FSAFM) 
in a magnetic field is a useful model system [3,4]: 

H = J J2 Sg-s$+J'Yl 4 -Sy-h^si. ( 1 ) 

(x,ff) [x,y\ x 



Here the S3 are spin-5" operators at the sites a: of a square lattice and the first 
and second sums are over nearest and diagonal neighbor pairs, respectively. 
This model is realized with S — ^ in Li 2 VOSi04 and Li 2 VOGe04 [5]. 
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2 We normalize the magnetization to \{M)\ < 1. 
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For J' — ^ the Hamiltonian (1) can be rewritten (up to a constant) as 

^4e{j3-k} (2) 

" a 

with L a = J2xea &x the total spin of a (four-spin) plaquette a and n — 2. In 
two dimensions the Kagome and checkerboard lattice and in three dimensions 
the pyrochlore and garnet lattices share this special form of the Hamiltonian. 

Consider first the classical variant of the model (2) with unit spins. Then, 
any state with L z a — jj, L x a = L y a = for all a is a groundstate, result- 
ing in huge degeneracies. Special collinear groundstates exist at special val- 
ues of the magnetic field (an up-up- up-down 'uuud' state for (M) — \ in 
the FSAFM, checkerboard and pyrochlore lattice and an up-up-down 'uud' 
state for (M) = | in the Kagome and garnet lattice). These collinear states 
have a particularly large number of zero modes leading to an exceptionally 
large number of low-energy excitations above them. At finite temperature, 
the probability function is therefore peaked close to the collinear states along 
the groundstate manifold leading to a selection of this submanifold by ther- 
mal fluctuations. The subsequent lifting of degeneracy goes differently for the 
FSAFM and Kagome AFM. In the former case a unique classical state is fur- 
ther selected [3], whereas in the latter case there is no such selection and a 
collinear spin-liquid state appears on the (M) = | plateau [6]. 

Semiclassical analysis [3] suggests that collinear states are also stabilized by 
quantum fluctuations, leading to a plateau in the magnetization curve of the 
quantum spin model. In fact, plateaux with (M) = | are expected in the 
quantum FSAFM, checkerboard and pyrochlore lattice and with (M) — | in 
the Kagome and garnet lattice even if they may not correspond to collinearly 
ordered states. An (M) = | plateau is indeed observed in exact diagonalization 
studies of the S — | Kagome lattice [7]. 

To test these general predictions, we have first performed Monte Carlo simu- 
lations of the classical FSAFM (see Fig. 3 of [3] for the magnetization curve 
at J' = j). Here, we present results for the static structure factor S a/3 (q) = 

y^J2x e ^ x (^g^f) i s the total number of spins) on a 24 x 24 cluster in 
Fig. 1. The top panel shows that at h — 4J non-vanishing harmonics in the 
longitudinal structure factor both at q — (ir, 0) and q = (n, ir) exist only for 
0.494 < y < °- 508 - Neither, a lower field of h = 2J, nor a higher field of 
h = 6 J induces a nonzero value of S zz (q) at these points in the Brillouin zone. 
The lower panel of Fig. 1 presents the transverse structure factor. In the re- 
gion of weak (strong) diagonal coupling J' one identifies Neel (striped) types 
of spin configuration. At h — 4J the corresponding peaks jump to zero at 
exactly the same points where the peaks in the longitudinal structure factor 
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Fig. 1. Longitudinal (top panel) and transverse (bottom panel) components of the 
static structure factor S a(S (q) for the classical FSAFM at T = 0.1J. 

appear, identifying the spin configuration in the intermediate region as the 
'uuud' state. For h = 2 J and h = 6J, no significant peak is observed in any of 
the structure factors for J' « ^, indicating disordered states above and below 
the uuud state. 

Finally, we study the S — \ FSAFM using Lanczos diagonalization. A detailed 
discussion [3,4] of the magnetization curves and the static structure factors 
shows that an (M) = | plateau exists in the region 0-5 < J - < 0.65 and that 
it corresponds to a state with uuud order. The finite-size analysis of states 
with (M) 7^ | suffers from the problem that data for different system sizes is 
available only for very restricted values of (M). Here we propose to interpolate 
the data for a given cluster size. The contour-lines which one obtains in the ~ 
(M)-plane are shown in Fig. 2 for 6x6 and 8x6 clusters. One indeed observes 
pronounced peaks on the 6x6 cluster at (M) = | which characterize the uuud 
state. The comparison with the data for the 8x6 cluster (which is available 
only for (M) > g « 0.70833) suggests that the peak in S zz (n,0) survives 
up to magnetizations as large as (M) 0.75 whereas S zz (7[,ir) presumably 
vanishes in this range. Although it remains to be clarified whether this (partial) 
order for (M) > | in the longitudinal components coexists with order in the 
transverse components (e.g. in S xx (n, it)), this signals a supersolid phase above 
the uuud state. This scenario is in sharp contrast to the classical FSAFM, but 
similar to the 2D S — | Shastry-Sutherland model where supersolid states are 
presumably realized above the (M) = | and (M) = | plateaux [2]. 

To summarize, thermal as well as quantum fluctuations stabilize a collinear 
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Fig. 2. Contour plot of S zz (ir,0) (top panel) and S zz (tt,tt) (bottom panel) for the 
S = \ FSAFM at T = 0. 

uuud state in the FSAFM for J' » f at (M) = \. For S = f, this leads 
to a pronounced plateau in the magnetization curve. In contrast to these 
similarities at (M) = |, the classical and quantum model presumably differ 
for magnetic fields above this uuud state: Disorder in the classical system 
versus a partially ordered, supersolid state in the quantum system. 
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